Optimum control of broadband noise by barriers based on sonic crystals 
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It is demonstrated that sonic crystals (periodic structures of sound scatterers) can be used to 
design acoustic barriers that attenuate efficiently broadband noise. Traffic noise is chosen here as 
an example in which our design procedure is applied. The structures consist of cylindrical units 
containing rubber crumb, a sound absorbing material. An optimization algorithm is developed to 
obtain the material distribution and the dimensions of the sonic crystal giving the best attenuation 
properties for this noise. The good agreement found between predictions and measurements for a 
barrier (3m height) characterized in a transmission room gives strong support to our proposal. 

PACS numbers: 43.50.Gf, 43.20.Fn 
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I. INTRODUCTION 

Attenuation of broadband noise is a topic of increasing 
interest because of the damage producing on human be- 
havior and social environment. Broadband noise is gen- 
erated by industrial machinery and many other products 
that we use in our modern lifestyle. Particularly, traffic 
noise is perhaps a paradigmatic example of broadband 
noise where different approaches for its attenuation have 
been extensively researched; see, for example the reviews 
in [l|, and references therein. Since the late nineties, 
two-dimensional sonic crystals were proposed [3-7] as an 
alternative to conventional noise barriers. Sonic crys- 
tals consist of periodic arrangements of scatterers that 
inhibit sound transmission for certain ranges of frequen- 
cies called bandgaps, as photonic crystals do with light. 
Optimization algorithms have been applied to this type 
of structures in order to enhance their acoustic perfor- 
mance. For example, Hakansson and coworkers applied 
genetic algorithms to design flat acoustic lenses @, de- 
multiplexers Q and many others acoustic devices. The 
same optimization procedure was recently used to en- 
hance the attenuation properties by means of the creation 
of vacancies in the periodic structures (T^]. However op- 
timization procedures recently studied in noise barriers 
design based on periodic structures have two main draw- 
backs. First, the barriers are very thick structures (i.e., 
consist of many layers of scatterers) making them prac- 
tically unfeasible for their construction due to their cost 
and space requirements. Second, the noise barriers that 
are placed at opposite sides of the traffic road are not 
taken into account. Under this circumstance multiply 
reflected sound between parallel barriers can cause a sig- 
nificant increase in noise in the screened area [111 ]. 

This work introduces a design procedure of sonic 
crystal-based barriers intended to attenuate efficiently 
broadband noise. The procedure is specifically applica- 
ble to traffic noise, which is a class of broadband noise 
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that has been very well characterized and widely studied. 
The barriers are based on the combination of two attenu- 
ation mechanisms; one is sound attenuation at the Bragg 
frequencies (due to reflectance) and the second is sound 
absorption by the material employed in the sonic crystal 
building-units. The barriers here introduced are made of 
only three rows of cylindrical units that contain rubber 
crumb, a sound absorbing material obtained from recy- 
cling car tires. An optimization algorithm is employed 
to obtain the optimum dimensions and material distri- 
bution in the scattering units of each row, as well as the 
distances between scatterers. A selected barrier has been 
constructed and experimentally characterized to support 
our proposal. Let us point out that the approach here 
applied to attenuate broadband noise generated by traffic 
noise can be easily extended to other types of broadband 
noise provided that its spectral profile is previously char- 
acterized. 

The paper is organized as follows. Section [TT| briefly 
describes the multiple reflection model applicable to the 
double-barrier configuration. Section [iTTl presents the op- 
timization goal function and the optimization algorithm. 
Results predicted by the optimized design are discussed 
in Sec. IIV1 The experimental validation of the algorithm 
is reported in Sec. |V] where we present a physical real- 
ization of a barrier of 3 meter height. Finally, Sect. IVII 
summarizes the work and gives conclusions. 



II. MULTIPLE REFLECTIONS MODEL 

Consider two barriers placed at both sides of a road 
where the traffic noise is generated. By assuming an in- 
coherent source, the problem of multiple reflections in 
the double barrier configuration can be approached with 
a ray model where the total energy at each point corre- 
sponds to the sum of the energy of every arriving path 
[12| . The simplified one-dimensional model is described 
in Fig. [TJ where the noise source (the car) radiates a 
plane wave with energy Eq . Each time the wave impinges 
a barrier, a fraction T of the incident energy is transmit- 
ted through the barrier, another fraction R is reflected 
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FIG. 1: Scheme of the multiple reflection model. The power in 
R defines the number of reflections experienced by the sound 
wave with initial energy Eq. 



and the rest 1 — R—T is absorbed. The effective transmis- 
sion coefficient, T e ff, is defined as the fraction of energy 
transmitted through both barriers: 
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This simplified model can be enhanced by considering 
the case of an infinitely long and straight road where con- 
tinuous traffic approximately behaves as a line source in 
such a manner that, at enough distance from the individ- 
ual sources, sound levels fall 3dB per doubling of distance 
as cylindrical wavefront does [l3j]. In other words, an at- 
tenuation proportional to 1/r (being r the traveled dis- 
tance) is finally obtained due to geometrical spreading. 
The corresponding transmission coefficient is 
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where £ is the separation between barriers. 

Another model is to consider road vehicles acting like 
individual point sources with spherical spreading; i.e. the 
sound energy is decreased by a factor 1/r 2 the coefficient 
transmission becomes: 
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In fact there are numerous traffic noise models, some of 
them considering the noise source as incoherent punctual 
sources and other as a line source [lij]. Because of this, 
the three canonical cases described by Eqs. |T] - [3] are 
analyzed here. 



III. OPTIMIZATION ALGORITHM 

We chose to maximize the so called insulation index 
for airborne sound (DLr), which is used in Europe as a 



criterion to classify the acoustical performance of traffic 
noise barriers. It is defined in the European Normative 
EN 1793-2 [H as 
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where Li is the normalized traffic spectrum, Ri is the 
sound transmission loss and i is an index indicating the 
eighteen standard third octave bands with frequencies 
from 100Hz to 5kHz. 

Parameters Li in Eq. Q take into account the human 
hearing response and emphasizes the frequencies where 
the traffic noise is more undesirable, having a maximum 
weight into the 1kHz band. The measurement of Ri is 
standardized in ISO140-3 [l6[ and it is briefly explained 
in Sect. El Note that EN 1793-2 is intended to charac- 
terize the intrinsic properties of the barrier, disregarding 
in situ conditions. Then, no discussion about diffraction 
by the top end of the cylinders or ground effects will be 
made in this work. 

The proposed structure to be optimized consists of a 
sonic crystal made of three infinite rows of cylinders, 
where each cylinder has an inner rigid core of radius 
(k = 1, 2, 3) and a layer of rubber crumb between the core 
and an external radius (see Fig. [2]) . Each row has iden- 
tical cylinders separated at a distance D and the three 
rows are separated by distances d\ and d,2, respectively. 
Therefore the optimization model involves 9 independent 
parameters. Note that cylinders are aligned in a square 
lattice, thus lowering the flow resistance of the barrier 
and making it partially transparent to light. For compar- 
ison purpose, we also studied a barrier in which cylinders 
at the middle row are displaced a distance D/2, being 
placed in the dashed circles in Fig. [5] Since this struc- 
ture forbids the light passing through, it will be called as 
opaque barrier. 

Calculations of barrier parameters are performed by 
applying the multiple scattering method. For computer 
performance reasons we use a 2D model, that is, infinitely 
long cylinders are considered. For technical details of the 
method the reader is referred to Ref. [l3 and references 
therein. Within this framework it is possible to calculate 
the reflectance Rsc and transmittance Tsc °f infinite 
sonic crystals by integrating the energy flux in a unity 
cell [3, resulting: 
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where u) is the angular frequency and 9q the angle of the 
impinging plane wave with respect to normal incidence. 
C~ and are the amplitudes of the reflected and trans- 
mitted v- modes which propagate at an angle V , obtained 
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FIG. 2: Parameters employed in the optimization of a noise 
barrier based on a sonic crystal consisting of three rows of 
cylindrical scatterers containing rubber crumb. Shadowed re- 
gions define the location of rubber crumb. The dashed circles 
indicate the positions of cylinders for the opaque (hexagonal 
symmetry) configuration. 



by the expression 
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k being the wavenumber and v an integer number. Ab- 
sorption by the sonic crystal barrier can be calculated by 
assuming energy conservation, so Age = 1 — Tsc — Rsc- 
The acoustic model used in describing the rubber crumb 
is based on the complex dynamical mass density and 
bulk modulus proposed by Johnson [l9j] and Stinson [2(| ■ 
A detailed theoretical and experimental work about the 
application of this model within the multiple scattering 
method has been recently reported in Ref. [2lj |. 

Multiple scattering method works with coherent waves 
with well defined wave fronts, making difficult to model 
the diffuse sound field required by the ISO140 norm [l6| . 
Since this type of field is the basis of several standard- 
ized measurements, analytical predictions have been de- 
veloped in control noise research. A conventional calcu- 
lation is based on the assumption that the angular distri- 
bution of incident sound field on the surface of the sample 
becomes fully uniform [22J . Consider a surface element of 
the sample dS and an isotropic incident sound intensity 
I. The fraction of acoustic power that arrives to dS from 
the direction defined by the solid angle dfl is 



Similarly, the total transmitted power is obtained as 
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being T(0) the angle dependent transmission coefficient, 
that is supposed to be independent of ip. The angle av- 
eraged transmission coefficient is given by the ratio 
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However uniform distribution does not fully reflect the 
actual sound field, so an angle dependency correction 
based on Gaussian distribution is applied [Ij, [24| . Hence 
the averaged transmission coefficient of the barrier is cal- 
culated as 



T av M = 
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and the same procedure can be applied to reflection and 
absorption coefficients, giving R av (u>) and A av (ui) from 
Rsc (w, 0) and Age (<*>, 8), respectively. 

Coefficients Ri in Eq. ((4]) are calculated as 1 /T x where 
T x is the direct T = T av or effective T eff , T e ' // , Tj? ff 
transmission coefficients obtained by applying T av and 
R av to Eqs. [IE]. These parameters are calculated at 
several frequencies in each third octave band and then 
they are integrated in order to get a single Ri per band. 

Due to the difficulty of differentiating the 9- 
dimensional objective function, we have employed the 
Nelder-Mead optimization method 25j which is based on 
the simplex algorithm. With the purpose to globalize the 
search, several initializations have been programmed in 
order to ensure that not local maximums are reached. 
Physical constrains such us positive dimensions or no 
overlapped cylinders need to be taken into account in 
the algorithm. Also practical constrains like maximum 
barrier width of 1 meter have been imposed. Maximum 
external radii of 10 cm and distance D at least four times 
the maximum radius are set in order to ensure a partial 
transparency of the barrier based on square symmetry. 
For the case of opaque barriers (hexagonal symmetry) 
the minimal distance between cylinders is set to 1 cm. 



IV. RESULTS AND DISCUSSION 



dPi = I cos(6)dSdfl = IdS cos(0) sin(0)dipd0, (8) 

where and tp corresponds to the azimuth and the incli- 
nation angles in spherical coordinates, respectively. The 
total arriving power is calculated by integrating Eq. [5] 
over the surface of a semi sphere 

/■7T /"27T 

Pi = / / IdS cos(0) sm(0)dipd0 = irIdS, (9) 
Jo Jo 



Barrier parameters resulting from the optimization 
process are shown in Table [H where a distance I =10m 
between parallel barriers has been assumed in Eqs. © 
and ©. ft is noticed that, for the case of opaque barriers 
(A symbols in Tabic HJ, the algorithm always converged 
to the maximum of data constrains (maximum external 
radius and minimum distance between cylinders); that 
is, trying to make the sonic crystal more compact. The 
same applies to the case of semi-transparent barriers (□ 
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TABLE I: Barrier parameters (see Fig. [2]) obtained from the 
optimization algorithm. Length dimensions are in cm. Last 
row contains the airborne insulation index DLr. Note that 
the highest quality barriers, class B3, according to the Euro- 
pean normative is achieved when DLr > 24rf_B[15l]. 
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FIG. 4: (Color online) Photograph taken inside the transmis- 
sion chamber. The barrier is schematically described in Fig. [2] 
and has parameters: r± — V2 = -R3 =10cm, m = Ri2 = 0cm, 
r»3 =4. 5cm, di =30cm, afe =50cm and D =42cm. 
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FIG. 3: Angular-averaged reflection T a „(w) and absorption 
A av {ui) of the opaque barriers optimized through the T and 
T e ff models. 



symbols) where minimum distance D and maximum ra- 
dius are also obtained. In addition distances d\ and di 
are practically the same for each transmission model in 
this topology. 

Looking closer to the parameters obtained, especially 
the internal radii, it is observed that cases based on ef- 
fective transmission T' e ^ and T"** tend to the values ob- 
tained for the case of simple transmission, T. This effect 
is due to the geometrical spreading of waves. This fact 
is easy to understand since multiple reflections are negli- 
gible when attenuation by wave propagation is large, be- 
coming important only the first incidence of the wave in 
the barrier. Note that differences in DLr between the T 
and T' e ff, T"^ ^ cases are mainly due to the attenuation in 
the first travel from the midpoint between barriers to the 
first barrier. Therefore, models T and T e ff can be seen 
as two extreme cases between which practical cases will 
stay depending on their attenuation by propagation and 
the distance between barriers. While in wide highways 



T' e fj and T"^ models will give parameters approaching 
those given by T, situations such as narrow roads or rail- 
way lines will make these parameters closer to the T e ff 
model. 



It is worth to note that for the case of simple transmis- 
sion T there are two insulation mechanisms available: re- 
flectance by the barrier periodic structure and absorption 
by the material of its building units. However the T e ff 
model only employs absorption since reflected waves al- 
ways reach the opposite barrier without being attenuated 
during propagation. The practical constrain of minimum 
distance D imposed to the optimization of the transpar- 
ent barriers has avoided the chance of using reflection as 
the principal attenuation mechanism (band gaps can not 
be formed on an efficient way). Therefore the optimiza- 
tion process has always chosen improving the absorption, 
thus resulting in barrier structures with similar parame- 
ters in the four transmission cases. On the other hand for 
the opaque barrier, where distances between cylinders are 
not so constrained, larger differences are found between 
the optimized parameters from the different transmission 
models. Note that in this configuration the barrier ob- 
tained from the T model needs larger rigid cores, espe- 
cially in the first row, while the barrier derived from the 
T e ff model requires more rubber crumb. This fact high- 
lights the differences between the transmission models in 
such a manner that in the first case a reflective barrier is 
obtained while the second corresponds to a barrier based 
on absorption phenomenon. The angular-averaged re- 
flection T av (uj) and absorption A av (oj) of these two cases 
are represented in Fig. [3j which shows how the barrier 
obtained from model T has a high reflective spectrum 
but a low absorption spectrum, while that resulting from 
model T e ff has an opposite behavior. Note that 
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V. EXPERIMENTAL CHARACTERIZATION OF 
AN ACTUAL BARRIER 

To validate our theoretical predictions, a 3m high by 
5m long barrier has been constructed and experimentally 
characterized in a transmission chamber (see Fig. U). 
Due to practical limitations, the optimized dimensions 
in Table U have been slightly modified, so that building 
parameters are R\ = R% = R3 = 10cm, Rn — Ri2 = 0, 
Ri3 — 4.5cto, d\ = 30cm, 6,2 = 50cm and D = 42cm. 
Theoretically the expected DLr of this barrier differs 
by less than 0.5 dB with respect to the optimum value 
reported in Table HI 

Coefficients Rj have been measured by following the 
1SO140-3 norm In brief, a loudspeaker is placed at 
the source room of the transmission chamber and gener- 
ates a white noise that becomes a diffuse sound field due 
to the multiple reflections on the room walls. Sound lev- 
els in source and receiver rooms (Li and L2, respectively) 
are acquired with moving microphones in such a manner 
that the sound field is averaged in time and space. Co- 
efficients Ri for each normalized third octave band are 
obtained as 



R(dB) = Lx- L 2 + Wlog 



(13) 



where S is the surface of the sample and A is the equiva- 
lent absorption area of the receiver room. The parameter 
A is obtained through the reverberation time T measured 
in the receiver room applied to Sabine's equation 



0.161/ 



(14) 



being V the volume of the receiver room. 

Figure [5] shows the Ri coefficients measured and cal- 
culated for models T (bold circles) and T e ff (bold trian- 
gles). According to the absorptive properties of rubber 
crumb, the responses of attenuation obtained increase as 
frequency does [2l|. As previously discussed, this barrier 
employs absorption as the main attenuation phenomenon 
while reflectance plays a minor role. Note that although a 
bandgap can be found in the reflectance spectrum, it does 
not appear if Fig. [5] because the angular and frequency 
average performed on the Ri coefficients smoothes the re- 
sponse. The agreement between theoretical predictions 
and experiment (hollow squares) is remarkable, being the 
obtained experimental value of DLr — 6.78dB, close to 
those predicted by the two models employed. Theory and 
experiment differ only at low frequencies where normal 
modes of the room become important, resulting in a not 
completely diffuse acoustic field. Moreover, as discussed 
above, both models produce similar curves due to the 
reduced reflectance of the barrier (see Eq. (JTJ). 

For comparison purpose, the case of a sonic crystal 
barrier made of only rigid cylinders (with the same ex- 
ternal radii) has been also considered and its coefficients 
are also depicted in Fig. [5J Note that its corresponding 
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FIG. 5: (Color online) Ri coefficients and quality fac- 
tors (measured and simulated) of barrier characterized in the 
transmission chamber. Coefficients for a barrier made of only 
rigid cylinders with the same external radii and for a flat bar- 
rier with the same quantity of rubber crumb per unit length 
are also depicted for comparison. 



quality factor has been strongly reduced in comparison 
with the rubber crumb barrier, where absorption is the 
mechanism leading to the broadband sound attenuation 
needed for traffic noise control. 

Figure[5]also shows the coefficients calculated for a flat 
panel of rubber crumb with thickness d = 20.9cm, thus 
having the same amount of this material per unit length 
than the measured barrier. These coefficients have been 
also obtained by using Eq. (1121) . in which the transmis- 
sion Tsc is replaced by that calculated for a slab made 
of a dissipative material (the rubber crumb) [26| 



T RC (0) 



-■yd cos((9) 
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(15) 



where 7 = a + jk is the exponential propagation in rub- 
ber crumb including attenuation effects and R the re- 
flectance between two semi-infinite layers of air and rub- 
ber crumb. Note that the DLr value of the flat barrier is 
slightly higher than that for the sonic crystal barrier, but 
this small difference is clearly compensated with the im- 
provement of flow resistance as well as aesthetic aspects 
of the sonic crystal barrier. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have used an optimization procedure 
to obtain optimal traffic noise barrier designs based on 
sonic crytals that attenuate efficiently broadband noise. 
The procedure uses the quality index DLr given by the 
European normative . Since standardized measures as- 
sumes an acoustic diffuse field, a theoretical framework 
to predict the behavior of sonic crystals made of rub- 
ber crumb in this kind of field is proposed. In order to 
consider realistic situations, multiple reflections between 
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parallel barriers have also been taken into account in the 
calculations, solving three canonical models of transmis- 
sion. It is shown that, in general, the optimization al- 
gorithm tends to make the barrier as dense as possible, 
although large differences can be found between results 
from each model of transmission. As a consequence, we 
conclude that schemes based on multiple reflections al- 
ways prefer barriers mainly absorptive. On the contrary, 
if multiple reflections are neglected, reflectance becomes 
a useful attenuation mechanism so reflective barriers are 
obtained. The experimental characterization of a pro- 
totype barrier constructed in a transmission room has 
given a strong support to the simulation algorithm here 
reported. Finally, our optimization algorithm can be ex- 



tended in designing sonic crystal-based barriers to atten- 
uate other types of broadband noise, as those coming 
from industrial machinery, with well defined (stationary) 
profiles. 
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